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Abstract 

It is experimentally observed that the width of the KNO multiplicity distribution 
— or the negative binomial parameter, 1/k — for pp collisions, in the energy region 
10 ^5 ^5 1800 GeV, is an increasing function of the energy. We argue that 
in models with parton or string saturation such trend will necessary change: at 
some energy the distribution will start to become narrower. In the framework of 
percolating strings, we have estimated the change to occur at an energy of the 
order of 5-10 TeV. 



1 The problem 



Experimentally, in pp and pp collisions, not only the average particle mul- 
tiplicity is increasing with energy but, as well, the normalised width of the 
multiplicity distribution is an increasing function of energy [1]. In other words, 
both (n) and D 2 /(n) 2 = ((n 2 ) — (n) 2 )/(n) 2 increase with energy. 

The increase of (n) is expected: more particles are produced when more energy 
is available. The growth of D 2 /(n) 2 is, however, not so obvious. The general 
argument seems to be that the increase of energy implies phase-space opening. 
There is more room for partons, strings, clusters or clans, more room for 
multiple collisions and different impact parameters contributions, more room 
for harder collisions and additional QCD branching [2]. In all these approaches, 
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if there are no limitations, fluctuations increasingly dominate over average 
multiplicity: D 2 / (n) 2 increases and, asymptotically, it may reach a constant 
value (KNO scaling [3]). This kind of behaviour, as pointed out in [4], has 
been observed long ago in cosmic ray physics [5], with the negative binomial 
parameter k stabilising at a value of the order of 2.5-3, for y/s > 1 TeV. 

However, if saturation phenomena occur [6,7], i.e., too many objects occupy a 
finite region of interaction, then fluctuations are severely limited and D 2 /(n) 2 
will start to decrease. On general grounds, this is not unreasonable. As the 
number of intermediate objects (or final particles) is bounded by y/s (conser- 
vation of energy), while the region of interaction is constrained by the stronger 
Froissart bound, In 2 y/s, there exists a real possibility for saturation. It is also 
well known that, at fixed energy, when triggering on rare events [8] or selecting 
central impact parameter, 6 = 0, collisions [9], the KNO distribution, due to 
the higher density, is narrower. 



2 Theoretical framework 

We shall work in the framework of the dual string model [10]. Saturation, 
in this case, results from clustering and percolation of strings. In a collision, 
N s strings are produced and they may form clusters with different number 
of strings (or different sizes). These clusters act as random, particle emitting, 
sources. Following [11], we write for the multiplicity distribution, P(n), 

P{n) = J W(X)P(n;Xn l )dX, (1) 

where W is a normalised weight function, X is a clustering variable, P{n; Xn\) 
a convolution function for a cluster of strings, fi\ being the single string aver- 
age multiplicity. We assume for P{n\ Xfi\) a Poisson distribution, this being 
consistent with the single string distribution of low energy e + e~ annihilations. 
Note that, as emphasised in [11], an equation similar to (1), with the same 
function W(X), can be written for the inclusive Pt distribution, f(px)- 

From (1) we obtain 
and 

1 

k 

Defining now N c as the number of clusters, N the number of strings in a 
cluster and N s the average number of strings, which must satisfy the sum rule 
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(N C )(N) = N s , we shall write 



X = (4) 

such that, finally, we obtain 

(n) = (N c )(N)n 1 =N s n u (5) 

and 

1 (iV 2 ) - (NT 

k = (N)> ■ (6) 

The relevant parameter in string clustering (and percolation) is the transverse 
density 77, 

I-®**, (7) 

where r$ is the transverse radius of the string (we shall use ro = 0.2 fm) and R 
is the effective radius of the interaction region, with R m 1 fm in pp collisions. 

It is the parameter r\ that controls clustering: when N strings cluster, the 
effective colour charge is not proportional to N but, due to the vectorial nature 
of colour summation [12], is proportional to F(r))N, where [13] 



1 - e _? ? 

F(V) = \j^ r , (8) 

such that, instead of (5), one should write 

(n) = F(rj)N s fii. (9) 

Equation (6) for 1/k remains, of course, unchanged. Note that F(r)) is a de- 
creasing function of rj. 



3 The toy model 



In order to show, clearly, the problem of the competition between the opening 
of new possibilities and limitations due to saturation, let us introduce a very 
simple, and known to everybody, toy model. Imagine that one has N s identical 
coins to distribute over M identical boxes. The reader may imagine the dis- 
tribution of N s partons or strings over the interaction region of size M, or R 2 
(in units of 7q). This is a purely combinatorial problem, without percolation. 

What is the distribution of coins in the boxes? In Ref. [14] it was shown that 
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Fig. 1. Plots of (a) the average number of strings per cluster, Eq. (11) and (b) of 
the inverse of its normalised variance, or k, Eq. (12), in the toy model, vs. the string 
density rj. 

in the thermodynamical limit only the transverse density 

" s j£- (10) 

is relevant, and that the first moments of the distribution P(N) of the number 
of strings per cluster (the probability of having a iV-cluster) 

(N) = ^- -, (11) 

V 1 l-exp(-?7)' v ; 



and 



K {N)2 3 (12) 

(N 2 )-(N) 2 1 - (1 + 77) exp(-?7) ' v ; 



If we use this toy model as input for W(X) in (1), and take, as mentioned 
before, the Poisson distribution for P(n;fii), then, see (6), one obtains 

K = k. (13) 



In Figure 1 we present plots of (11) and (12). The average multiplicity is, 
naturally, a growing function of the density. The parameter k goes to infinity at 
small 7] (one just has 1-clusters), decreases and reaches a minimum, k ~ 3.3, at 
77 ~ 1.8, then moves to infinity again as rj — > 00 (the distribution is dominated 
by iVs/M-clusters). We have checked by Monte Carlo simulation that the 
distribution P(N) of the model is similar to a gamma distribution. 
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4 The physics 



In order to make contact with real physics, we need to estimate the dependence 
of the parameter k, experimentally measured in pp collisions as a function of 
the energy, on the density parameter rj. In other words, we have to find the 
energy dependence of 77. As i], see (7), depends on N s and R 2 , we need to find 
the energy dependence of iV s and R 2 . 

For the pp average charged particle multiplicity in the dual string model, we 
shall write [15]: 



where fi\ is the single valence string average multiplicity (assumed constant), 
a the sea string multiplicity reduction factor (a < 1). Notice that in (14) one 
always has 2 valence strings and N s — 2 sea strings. At low energy, y/s pa 
10 GeV, experimentally (n) pa 4, and theoretically N s pa 2 (the two valence 
strings) and of course F(rj) ~ 1 and fi\ ~ 2. Notice that (14) is a natural 
generalisation of (9) when there are two kinds of strings. 

As the slope parameter B(yfs) is a direct measure of the square of the inter- 
action radius R 2 (y/s), we write, from (7), 



with c ~ 1.6 GeV 2 fm 2 , such that at y/s pa 10 GeV, it! pa 1 fm and B pa 
10 GeV" 2 . 

By using equations (14) and (15), and experimental information on (n(y/s)) 
[16] and on B(y/s) [17], we can in principle construct N s (y/s) and r)(y/s): from 
v{V~s) we obtain k(y/s), Eq. (12). 

In Fig. 2a we show our fit to B(y/s) and in Fig. 2b our fit to (n(y/s)). These 
parameterisations are then used in Eq. (15) and Eq. (14). 

Intuitively, one expects the number of strings, N s , to increase with \/s. The 
results from our fits are shown in Fig. 3: N s increases with energy and, within 
a factor 2, is not substantially different from Monte Carlo direct calculations 
[23]. Naively, as r] ~ N s , one also expects i] to grow with energy. However, 
as the interaction area also grows with energy, the situation is not so clear. 
From (14) and (15) it is easily seen that the 77 dependence on y/s, r/(y / i), is 
determined by the ratio (n( A /s))/i?(v / s). In Figure 3b we present r)(y/s); in 
the figure is shown as well the percolation threshold: r/ c ~ 1.15. 



(n(v^)> = F^n, 2 + (N 8 (yfi) - 2) a , 



(14) 




(15) 
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Fig. 2. (a) Slope parameter B (GeV 2 ) vs. cm. energy (data points from 
[17]), parametrised by B(y/s) = 9.52 + 0.671 In ^/i + 0.0424 In 2 ^/i- (b) Av- 
erage multiplicity (data points from [16]) vs. cm. energy, parametrised by 
(n(y/s)) = 5.99 - 2.46 In + 0.987 In 2 ^/i. 

Finally, in Fig. 3c, we show the dependence of k on energy, see (12) and Fig. 3b, 
in comparison with experimental points. If the "coins in boxes" model reflects 
the physical situation, we expect k to reach a minimum at ^/s pa 6 TeV, in the 
LHC energy region. The behaviour we predict is qualitatively different from 
what is expected on the basis of extrapolations of other phenomenological 
models: for instance, consider the Wroblewski regularity [19]: D = a(n) + b; 
this leads to an asymptotic ((n) — > oo) KNO scaling, thus to a levelling off of 
k (see dashed line in Fig. 3c.) Other proposals (e.g. k~ l = a + b\n(n) in [20], 
and the Pythia Monte Carlo [22], respectively dotted line and cross in Fig. 3c) 
give a monotonically decreasing k; for a discussion of the situation in which 
k continues to decrease at high energy and reaches values smaller than 1, see 
[21]. 



5 Conclusions 



We have argued that the energy behaviour of the parameter k (or (N) 2 / D 2 N ) 
of the multiplicity distribution in pp collisions is sensitive to the energy de- 
pendence of the parton/string transverse density rj. If impact parameter over- 
crowding never occurs (rj is always small) k will decrease with energy, reaching, 
eventually, some asymptotic constant value. If overcrowding takes place (rj in- 
creases rapidly with energy) then k necessarily will change its behaviour and 
start increasing with energy. The turning point may occur at an energy of the 
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Fig. 3. (a) N s vs cm. energy; (b) r\ vs cm. energy; (c) k vs cm. energy; experi- 
mental points from ISR and UA5 (filled points, [16]) and E735 (open squares, [18]). 
The values of the parameters used [see (14) and (15)] are: a = 0.45, fi\ = 2.9 
and c = 1.9 GeV 2 fm 2 . The dashed line is a fit using Wroblewski regularity, 
D = 0.59(n) - 1.4; the dotted line is a fit using k' 1 = 0.154 In (n) - 0.25 (both 
fits have been carried out considering only data with y/s > 30 GeV.) The cross is 
the prediction by the Pythia Monte Carlo with the most recent parameters optimi- 
sation [22]. 

order of 5-10 TeV. 

It was shown in [25,11] that the parameter k of the multiplicity distribution 
is the parameter k of the p^-distribution: (1 + bp^)~ k . This parameter k is 
also a decreasing function of energy in the y/s < 1.8 TeV region. In Figure 4 
we have compared experimental data on distributions to the mentioned 
parameterisation, with k given by the multiplicity distribution values at the 
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Fig. 4. px distributions (not normalised) at several cm. energies, (from [16,24]) 
fitted with the power-law A(l + bp\)~ k , using, for each energy, the respective value 
of k obtained from the multiplicity distribution (from Fig. 3), A and b being free 
parameters. 



same energy. The agreement is reasonable, and improving with energy. 
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